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Abstract 

We re-examine the possibility to measure the photon polarization in B —> K* 7 decays, via 
decays in which the photon subsequently undergoes nuclear conversion to a lepton pair. We obtain 
compact expressions for the full decay-plus-conversion amplitude. With these results we show that 
interference between the B —> (K* —y Kir)'y decay and the 7 N —> £ + £~N conversion permits 
both the ratio and relative weak phase between the left- and right-handed photon amplitudes to 
be probed by an angular observable, constructed from the final state dilepton, kaon and pion 
kinematic configuration. Exploiting this technique will be experimentally challenging. However, 
we present special kinematic cuts that enhance the statistical power of this technique by an 0(1) 
factor. We verify this effect and extract pertinent angular kinematic distributions with dedicated 
numerical simulations. 
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I. INTRODUCTION 


Measurement of the photon polarization in b —> sy radiative decays has been of long¬ 
standing interest. Within the Standard Model (SM), flavor-changing electroweak interac¬ 
tions maximally violate parity, so that one expects the fraction of left-handed photons in 
b — y sy processes to be order unity, up to small corrections arising from either the non¬ 
zero strange quark mass or from higher order QCD contributions. In contrast, certain New 
Physics (NP) scenarios may generate b —» sy operators of comparable size to the SM terms, 
but with exotic parity structure, significantly modifying the expected ratio of left- versus 
right-handed photons - the photon polarization ratio. Measurement of this ratio therefore 
has the potential to test the parity structure of b —> sy operators against SM expectations, 
as well as either constrain or detect the signatures of such NP scenarios. 

The b — )■ sy photon polarization ratio may be measured via various different approaches. 
Dominantly right-handed photon production in resonant B —>- (A\(1400) —>- Knn )y gener¬ 
ates an up-down asymmetry of the photon momentum with respect to the K\ —>■ Ktitt decay 
plane [1, 2]. This up-down asymmetry was recently measured by the LHCb collaboration [3]. 
However, a theoretical prediction for the asymmetry is not yet available, so that the photon 
polarization ratio cannot yet be extracted from these results or compared to SM expecta¬ 
tions. Along similar lines, the photon polarization ratio may also be probed by measuring 
the spin fraction of A’s in unpolarized A& —)■ Ay decays [4], or by measuring an angular 
asymmetry between the A& spin and the outgoing photon momentum for polarized Aj, —>■ Ay 
[5]. Other methods look for time dependent CP violation induced by mixing of B —> K* y 
and B —> K* y, which is proportional to the photon polarization ratio [6, 7]. Additionally 
one may probe the polarization ratio by looking for asymmetries in angular observables in 
resonant B —>■ Kir£ + £~ or B — > nnt + £~ [7, 8], or look for transverse asymmetries in the 
dilepton invariant mass for non-resonant B —>- Kti£ + £~ [9]. 

In this work we focus on the B —> ( K* — > K + n~ )y process, in which the on-shell pho¬ 
ton subsequently undergoes Bethe-Heitler (BH) nuclear conversion inside the detector to 
a lepton-antilcpton pair. The cross-section for BH pair conversion of ~ GeV photons is 
approximately two (eight) orders of magnitude larger than the Compton (Rayleigh) scatter¬ 
ing cross-section ([10]; see chapter 32), so that to an excellent approximation the emitted 
photon does not decohere before conversion. The B —> K* y photon polarization ratio, r, is 
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precisely defined via the amplitude ratio 

i(t+6) = ->• l) 

~ A(B -y K* 1r ) ’ 

in which 0 (5) is a relative weak (strong) phase. In the SM, r is expected to be at most 
~ A qc d /mb [11], while the weak phase 0 is suppressed by A qc d /mb V( ife V(* s /V) 6 V)* -C 1. 
Hence measurement of not only r, as discussed above, but also 0 may test SM expectations: 
Measuring either r or 0 ~ 0(1) would be highly suggestive of NP effects. 

Measurement of r via BH conversion was first considered in Ref. [12]. In that analysis, 
the B meson was assumed to be at rest relative to the conversion nucleus. Further, the 
conversion itself was assumed be to a perfect linear polarizer of the photon, so that the 
conversion leptons and photon are constrained to be coplanar. However, in practice the 
B meson is typically at least semi-relativistic, and typically an 0(1) fraction of conversion 
events have non-negligible acoplanarity (see e.g. Ref. [13]). This leads to a richer phase space 
structure for the outgoing conversion leptons, kaon and pion. Moreover, interference effects 
between the B —> K *7 decay and the BH conversion amplitudes were not included. The key 
motivation to reconsider the above analysis, then, is to include these R-boost, acoplanarity 
and interference effects. We exploit recent compact results for BH conversion spin-helicity 
amplitudes [13] to construct the full B —> ( K* —> K + TT~)('f —>• e + e _ ) amplitude, and show 
that interference between its decay and conversion components permits both r and (j) to be 
probed by the kinematic configuration of the final state conversion leptons, kaon and pion. 
In particular, we develop an angular observable that probes both of these parameters. This 
new analysis admits arbitrary boosts of the parent B meson relative to the BH conversion 
nucleus, and includes lepton-photon acoplanarity, which turns out to play an important role 
in enhancing the r- and ^-sensitive interference effects. 

Performing an experiment to measure r and 0 with this technique will be challenging. 
In the first instance, precise reconstruction of the leptonic momenta is required, which can 
only be achieved, even in principle, if the leptonic opening angles are larger than the angular 
resolution of the detector. Typically, the leptonic opening angle after BH conversion is 
&u ~ m e /E 1 rv_/ 10 4 , for typical photon energies in a B —> K *7 decay with a semirclativistic 
B. Specifically, for a photon of energy < 5 GeV, the probability for 9u > 10 -4 (10~ 3 ) is 
~ 98% (43%). Hence exquisite angular resolutions will be required. A further complicating 
factor is the multiple rescattering of the leptons in the detector material after conversion. 
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The rms rescattering angle in matter is ~ (13.6 MeV/i%) y/ x/xq [10], where x/xq is the 
path length inside the detector in units of radiation length. For x/xq ~ few % - a typical 
value - the rescattering angle is comparable to the typical opening angle, 9a- Finally, the 
probability of photon conversion itself is typically low at current and planned B factories, 
being at most of order a few percent. This probability depends mildly on the detector design. 
For example, it is approximately 3%, 2-3%, and 6 % at BaBar [14], LHCb [15], and Belle 
II [16] respectively. For all these reasons, this technique will likely only be feasible with a 
dedicated detector element that has a large scattering length, e.g. a gaseous TPC [17]. 

For these reasons, in this work we shall restrict ourselves to a thought experiment-type 
approach. That is, we develop explicit analytic expressions for amplitudes and observables 
with respect to the underlying B —> ( K* —> K + 7 r _ )(7 —> e + e~) process alone, but do not 
include smearing from leptonic rescattering and limited kinematic resolution, or realistic de¬ 
tector simulations. Based on the results of this work, future studies may perhaps incorporate 
these latter effects. 


II. AMPLITUDES 

A. Amplitude factorization 

Keeping operators up to dimension five, the effective theory of interest for B — > K* 7 and 
B — > K* 7 decays may be written in the general form 

As = JiiBtA'ye'- + g ± fit + g,B'K^ + g^K; w F “*, (2) 

where X tiu = e Ml/pt 7 X p(T /2, and the dimensionful couplings generically contain relative strong 
phases. We consider only K* decays to charged pseudoscalars, i.e. K* —> K + tt~ or K* —> 
K~ 7 r + . The sign of the pion or kaon charge therefore tags the K* versus the K*, and hence 
tags the parent meson as either a B or B, up to electroweak loop suppressed corrections. 
Hence we neglect interference effects from B-B mixing. 

We assume the conversion nucleus is spin-0, e.g. a 28 Si nucleus, which is the dominant 
silicon isotope. The external quantum numbers for the full B —> ( K* —Ji + 7 T - )(7 —>- e + e _ ) 
helieity amplitudes are then just the spins of the electron and positron, denoted r and 
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s = 1,2 respectively. It follows that the full helicity amplitudes 



nearby to the K* Breit-Wigner peak. Here A = ± (k = ±, 0) is the helicity of the photon 
(/!*), and (p M ) is the photon ( K*) momentum; P (P') denotes the incoming (outgoing) 
nuclear momentum, with nuclear mass P 2 = P 2 = M 2 ; p± (pK,n) denote the momenta of 
the leptons (kaon and pion); and finally q denotes the momentum exchange with the nucleus. 

Momentum and angular momentum conservation in the B —> K* 7 process ensure that 
As^K*-y must annihilate the longitudinal component of the K* propagator. Applying the 
polarization completeness relation for p 2 7 ^ 0 , 


e Z(P) eK Z(p) = -9^+PnPv/p 2 , 

/■c=d=,0 

we may factorize the full helicity amplitude into three helicity amplitude factors, 

i 


M r , = 


p 2 — m ? K * + im-K* T k* 


[Ab-> K*~/] XK [Ak*^K+it- ] h ' [^Bh] rs ■ 


\,K 


BH 

Similarly for the B —» (K* —> K~7T + )( r y —> e + e _ ) process 

i 


M rs = 


p 2 — m 2 K * + irriK*^ k* 


^ [A b ^r* 7 ] Ak [A** 


-kFr- 


] k [Abh] A s 


\,K 


( 4 ) 


(5a) 


(5b) 


Here and hereafter we neglect the mass splittings of the B-B and K*-K* systems, and denote 
the masses (momenta) of both CP conjugate states by m B and m K * ( p B and p) respectively. 


B. Kinematics 

The amplitude factors in eqs. (5) are Lorentz invariants, and are naturally expressed with 
respect to kinematic coordinates that are defined in different frames. That is, the full 2 —>■ 5 
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phase space - the full coherent process is BN —> K + n~£ + £~N, for nucleus N - is factored 
into two 1 —» 2 phase spaces, corresponding to B —y K* 7 and K* —y K + n~, and one 2 — y 3 
phase space for the BH conversion. I 11 general, a 2 —>• 5 phase space is fully specified by 
eleven coordinates. Two of these are the photon and K* invariant masses. The former is 
fixed to k 2 = 0 for an on-shell internal photon. With regard to the K* invariant mass, it is 
convenient to define hereafter the dimensionless quantity 

s = p 2 /m 2 K *. (6) 


In the narrow width limit, the Breit-Wigner factor 


1 

|(s - 1 )m 2 K , + im K *Y K * | 2 


1 

^-p— 

m K *i K * 


<5[(s - 1)774,] • 


( 7 ) 


That is, the narrow width limit corresponds to an on-shell K*. Since, however, the K* has 
a finite width - Yk*/ttik* ~ 5% ~ and need not be precisely on-shell, we shall treat s as a 
phase space variable: The Breit-Wigner ensures s is typically nearby to the K* mass shell 
up to the Ii* width, i.e. s ~ 1 up to variations ~ Y k */ijik*- 

In the case that the photon conversion material is cold, the lab frame coincides with the 
frame in which the BH conversion nucleus is at rest. The following choices, shown in Fig. 1, 
for the remaining nine coordinates then prove convenient for the construction of compact 
and intuitive results: the lepton polar and azimuthal angles 6± and <f>± and the energies E±, 
defined in the nuclear rest frame - the lab frame - with respect to the photon momentum 
and the K*- 7 decay plane, defined by p and fc; the photon polar angle, 0 7 , defined with 
respect to the nuclear momentum, P, in the B rest frame; the K polar and azimuthal angles 
6k and 4>k, defined in the K* rest frame with respect to the photon momentum, k, and the 
plane defined by P and k in that frame. Note that the K*- 7 decay plane is invariant under 
boosts between the lab, K* and B rest frames, and therefore equivalent to the plane defined 
by P and k in either the K* or B rest frames, as shown in Fig. 1. 


C. Helicity amplitude factors 

With these choices, we now proceed to explicitly compute the amplitude factors [Mb_^* 7 ] Ak , 
[M K *^. K + n -] K , and [Mbh]^- Applying a light-cone decomposition to p, we define its associ- 
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FIG. 1. Kinematic configuration and coordinate choices. B momentum is denoted by ps , and 
azimuthal angles are defined with respect to the K *-7 decay plane (blue). This plane contains p 
and k (P and k) in the lab frame (K* or B rest frames); momenta lying in this plane in each 
frame are shown in gray. Left: Lepton polar angles 9± and azimuthal angles 4>± in the lab frame. 
Middle: 9k and (j)x polar angles in the I\* rest frame. Right: The photon polar angle, 0 7 , in the 
B rest frame. 


ated null momentum with respect to the photon, i.e. 


p^=p^~ 


p 2 k fl 
2p ■ k ’ 


and make the polarization gauge choices 


( V) = ± 


(kA\a^\p T ) 
V2(k*\p±) ’ 



± (k^\a^) 

^(k^lkA) 


( 8 ) 


(9) 


for k an arbitrary null reference momentum. From the effective theory (2) one may then 
show that the B —» K *7 and B —y K *7 helicity amplitudes are 


[Ab^k*^ = (g\\ ± ig±)(m% - m 2 K ,s ), 
[Ab^k^ = (g\\ ± igx)(m% - m 2 K *s ), 


( 10 ) 


and [MB^K*-y] T± = [Mb^k*^ 0 = 0. Note that, by definition, A{B —» K*<^r,l) = 
[Mb_*.x* 7 ] ±± - It follows from eq. (1) and its CP conjugate that 

911 — — = re i{5+(t,) and ? + l ' 9± = . ( 11 ) 

fl'n + ig± g 11 - ig± 
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Hence 


[A B ^K*-y] ++ = (g\\ + ig±)(m% - m 2 K ,s ), 

[A B ^K*^r~ = re l{ - 5+<t>) (g\\ + ig±)(m 2 B - m 2 K ,s ), 

[Ab^k*j} ++ = re t{ - s ~^{g\\ - ig±)(m% - m 2 K ,s ), 

[A S ^k*t\— = i.g\\ ~ ig±)(m 2 B - m 2 K As) . (12) 


The reference gauge momentum k in eqs. (9) is so far arbitrary. However, a particularly 
convenient choice is 


W = 


2 P-k 

M 2 


P p - k» , 


(13) 


where P 2 = M 2 is the nuclear mass. In the lab frame - the nuclear rest frame - this 
choice (13) ensures that for k M = E^(l,k) then simply k 11 = E y (l, —k). We assume the 
nuclear scattering is coherent and quasi-elastic, i.e. that P' 2 = M 2 - equivalently q° = 
P'° — P° = q 2 /(2 M) in the lab frame - and that the outgoing nucleus is non-relativistic, so 
that the momentum exchange with the nucleus |qr| <C M (see Refs [18-20] for a review of BH 
conversion). With these assumptions and the choice of k in eq. (13), the BH spin-hclicity 
amplitudes collapse to a simple form in the limit that the polar angles &± <C 1 and Lorentz 
factors 7 ± = E±/m S> 1, where m is the lepton mass [13]. At leading order in these limits, 


[A. 


BH 


~ e 3 au 


(14) 


with 


a 


li 


= -(«+)* = 2^2^Z 




q 2 V 1 + 7+0+ 1 + l-0 2 -J ’ 


- _ ,f / o„ V^O? 2 ) 7+ ( 7+0+e i<t>+ , 7 -8-e~ 

— +(«2i) — i 2 y 27 _|_ 7 _--;- I “7—;-ttwi b 


q 2 7 + + 7_ V 1 + 7+0+ 1 + 1-9- J ’ 


a 22 ~ ( a n) — 0 1 


in which 


(15) 


9 

\ m 

r i 

11 

+ TT 

— + — 

/ 4 

L7+ 

7-J 


(16) 


Here G(q 2 ) is the BH quasi-elastic form factor for the photo-nuclear vertex [19, 20], 


G(q 2 ) — M 2 a 4 q 4 /(1 — a 2 q 2 ) 2 , 


( 17 ) 















in which a = 184.15(2.718) _1 / 2 Z _1 / 3 /m and Z is the atomic number of the nucleus. This 
form factor encodes electronic screening of the nucleus, and regulates the 1/g 2 pole in the 
amplitudes. From eqs. (15), one sees that the BH amplitudes are maximal at 7 ±9± ~ 1. That 
is, the typical lepton polar angle 9± ~ m/E±. It also follows that the typical momentum 
exchange — q 2 ~ m 2 , i.e. |qr| 2 M 2 in concordance with our assumption of non-relativistic 
scattering. 

It now remains to compute the K* —> K + tt~ and K* —> K~tt + helicity amplitudes. Only 
transverse K* modes are generated by the B —> K* 7 amplitude. Since e^-* • p = 0 and 
P = Pk + P-k, these amplitudes must therefore take the form 


\A k .^ k +v-Y = - K) . 


(18) 


where gx* is a dimensionless coupling, containing a strong phase. There are no other weak 
phases as K* decays strongly to Ki r. Under the polarization conventions (9), and computing 
in the Ii* rest frame defined by Fig. 1, the helicity amplitudes are just spherical harmonics 


0 ±i<j>K 


[Ak*^k+it-} (^k,<Pk) — PkkPktt sin 6k , 


with the momentum 


PKn — 


[m 2 K ,s - ( m K + m n ) 2 ] 1/2 [m 2 K «s - ( m K - m,) 2 ] 


S l Affix* 


\ 211/2 


(19) 


( 20 ) 


Under CP, note that the amplitude transforms as 


[CP Ax*^K+Tr-} ± (9x, 4 >k) — [A-K*^ K - 7r+ ] T (d K , 4 >k) — [AK*^K+n-] ± (^K, ~4 >k) • (21) 


That is, defining (f)x and 9k with respect to K* rest frame just as for the K* in Fig. 1, then 


0 ±i(f>K 


M]T^k-*+] (9k,4>k) = -^-<Jk*Pk* sin Ok ■ 


( 22 ) 


D. Full Amplitude 


Applying all the results (12), (15), (19) and (22) to eqs. (5), and defining |a 
^ ~2 Xrs | cx^ s | 2 , the unpolarized square amplitudes 


\M \ 2 

\M \ 2 


A(r ) sin 2 0 K 
A(r ) sin 2 9 K 


m + 


a + 


8 r 


1 + r 2 
8 r 

1 + r 2 


r Re 


r Re 


«r 2 « 21 e i(0+5 - 2 ^ ) 


ar 2 «2-ie i( ^- 2 ^ } 


(23) 
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in which 


A(r) = j( 1 + r 2 ) 

A(r) = j(l + r 2 ) 

Eqs. (23) compactly express the unpolarized square amplitude for the full B —* ( K* —> 
Jl + 7 T _ )(7 —> e + e~) process in terms of just the BH conversion helicity amplitudes (15) and 
trigonometric (exponential) functions of the kinematic observables 6k (4>i<)- The dependence 
on the parameters r, (j) and <5 is explicit and elementary. 

As a cross-check of these results, we provide an alternative and more traditional deriva¬ 
tion of the square amplitude in Appendix A, via construction of linearly polarized pho¬ 
ton BH amplitudes. The consequent result (A4) and the square amplitude (23) are in 
excellent numerical agreement in the 7 ± 1 and 9± <C 1 regime, applicable to the 

B —> ( K* —>• A' + 7 t _ )(7 —>• e + e _ ) process. Note that the result (A4) does not incorporate 
these approximations, so that the compact and explicit eqs. (23) are strictly an approxima¬ 
tion to eq. (A4). 


Iflji +ig±\ 2 \gKn\ 


2/2 2 -\2 
2PKn( m B- m K* S ) 


| m 2 K ,(s - 1 ) + im K *Y k 

\gw-ig±\-\gKn\ 2 2 , 2 2 . 

I 2 -F- \2PKn( m B - m K*S 

\ m K* I s — 1 ) + ^K* 1 K* | 


(24) 


III. OBSERVABLES 

A. Differential rate 

Making use of the explicit r and (j) dependence in the square amplitude results (23), we 
may now proceed to extract r and (j) sensitive observables. First, however, we construct the 
full differential rate. The factorization (5) ensures that the phase space with an on-shell 
internal photon may be partitioned into a B —> (K* —y Ktt )y cascade decay and a 7 N —y 

DU 

i + £~N conversion. That is, the differential rate for the full B —> ( K* —> K + TT~)('y —> e + e ~) 
process can be written as 

dn=\M\ 2 dV B dV B u (25) 

where dV B (dV B h) is phase space factor of the decay (conversion). Note that dlZ here has 
the dimensions of a cross-section times a partial width. 

Each phase space factor is Lorentz invariant, and are naturally computed in different 
frames, as shown in Fig. 1. Computing in the B rest frame followed by the K* rest frame, 
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the phase space factor for the cascade decay is 


cIVb — 


1 d 3 p 7r d 3 p K d 3 k 4 


27tib ( 27 r ) 5 2 E n 2 Ek 2 E, 


5 (pb — k — p K — p „), 


7 


m K *p Kn m 2 B -m 2 K Es 

—y ittt ;- - - _ ———— d\ IrCI COS u^CLS , 


( 26 ) 


s 1 / 2 32(27r) 4 m^ 

performing all trivial integrals, including over the overall azimuthal orientation of the K*- 
7 -N plane. Similarly, for the Bff conversion, computing in the lab frame, 


dV RH — 


1 d 3 p + d 3 p_d 3 P' 4 , 

o (P + k — P — p + — p _), 


->■ 


2M2E 1 (2vr ) 5 25+ 25_ 25' 
^ 2 -A 2 


_64(2tt) 5 M 2 ^ 


dQ.dQ_dA. 


(27) 


ffere the lepton momenta been approximated in the measure via yjE± — m 2 ~ E±. We 
have further defined 


£ = 


E+ + E_ 
2 



(28) 


and enforced non-relativistic nuclear scattering, which implies E' ~ M or equivalently E 1 ~ 
E + + E-, up to q° = q 2 /2M m corrections, ffence to an excellent approximation £ is 
half the photon energy in the lab frame. Moreover, note that m < E± < E 1 implies that 


A G (m — £, £ — m). 


(29) 


At an e + e“ 5-factory, such as Belle or BaBar, the e + e _ —> T —>■ BB production factorizes 
from the subsequent B decays, because the B is a pseudoscalar. In this type of collider, 
the rapidity of parent B meson in the lab frame has a known prior probability distribution, 
fs(p)dp, determined by the collider configuration, and enters as an independent phase space 
factor in dlZ. In Appendix B we include a derivation of the B rapidity pdf (B 6 ), for an e + e~ 
machine. We shall restrict ourselves hereafter to the case that fs{p) is known. In this case, 
note that the energy £ is fully specified by p, s and 0 7 , viz. 

2 9 

jYl _ 777 / S 

£ (p, dry, s) = ^- - —- ( COSll P + COS dry sinh p ) , (30) 

so that the lepton energies can be expressed in terms of p, s , and A, via eqs. (28). In 
our discussion of the kinematics above, it was convenient to express the amplitudes in terms 
of the ten phase space coordinates s, d 7 , E±, kl K and £l±. We see here, however, that for 
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the differential rate, it is more natural to choose rj and A as phase space coordinates rather 
than E±. Combining the above results together, the full differential rate 


dll 


1 Ml 2 fM) mK * PK7T 

m 2 B — m 2 K *s 

~£ 2 (v,^,s) - A 2 ’ 

r u l l B yl) - 1/2 

4(47r) 9 m^ 

M 2 £(rj, dry, s ) 


di),Kdfl+dQ-dAd cos O^dsdr). 

(31) 


B. Polarization and weak phase observables 


Let us now define a further change of azimuthal angular coordinates, modulo 2n 

0 = 0+ + 0_ + 2 (p K , 

0 = 0+ + 0_ - 2 <p K , (32) 

= <t>+ ~ <t>- ■ 

The angle p encodes the acoplanarity of the leptons with respect to the photon, with copla¬ 
narity corresponding to p — n. Note that <pK and (p± are defined with opposite orientations 
around the photon momentum direction k (see Fig. 1). For coplanar conversion leptons 
and a stationary B in the lab frame, ip then corresponds to the relative twist between the 
positron-electron conversion plane and the K-tt decay plane. Similarly, ip would then cor¬ 
respond to the averaged orientation of the positron-electron conversion plane and the K-tt 
decay plane with respect to the K *-7 decay plane. 

From eq. (16), the momentum exchange has the form q 2 oc 1 + ( cos p, with ( < 1. It 
follows from eqs. (15), (23) and (31) that the differential rate can be written in the form 


dll = ^ afc cos k (p) S J Ai + A 2 cos(ip + p — 0 — 5) + 

k 

A 3 cos(0 — p — 0 — S) + A 4 cos(ip — (p — 5) , (33) 

where a*, and Ai are purely functions of the phase space orthogonal to 0, 0 and p. That 
is, 0 ^ and Ai are functions of s, ip d l K ,± and A. Integrating over all phase space except 
dip, we see that the marginal differential rates for B —)■ ( K* —> K + vr _ )(7 —> e + e _ ) and 
B —>■ (K* —y K~7T + )( r y —> e + e ~) must respectively have the form 


dll 

dip 

dn 

dip 


n 
20 L 
n 
20 


2 r 


1 + r 2 
2 r 

1 + r 2 


C cos (0 — cp — h) 
C cos (ip — cp + S) 


(34) 


12 

















Eqs. (23) tell us that the cosine coefficient C arises from a ratio of BH interference terms 
to the BH squared amplitude, and is therefore independent of r or </>. In other words, this 
coefficient is the same for both of the CP conjugate processes, and is B —>• K *7 operator 
independent. We have chosen the relative sign of C in eqs (34) to anticipate the choice 
that ensures C > 0. Further, we have chosen the normalization of C in eqs (34) to ensure, 
via positive semi-definiteness of dlZ/di/), that \C\ < 1, noting that the r dependent factor 
2 r / ( 1 +r 2 ) < 1 for any r. The coefficient C may then be interpreted as the maximum possible 
ratio of the amplitude of dlZ/d 1 / oscillations to their average value, lZ/2n. Hereafter we call 
this ratio the relative oscillation amplitude. 

Eqs (34) are the main results of this paper. Once the coefficient C is computed, then 
measurement of the relative oscillation amplitude in dlZ/di/> permits extraction of r up to 
the two-fold ambiguity r +• 1/r. Further, measurement of the average phase shift (phase 
shift difference) between dlZ/di/ and dlZ/di/ permits extraction of the weak (strong) phase (j) 
(h). Equivalently, one may construct two forward-backward type asymmetries. Defining the 
four quadrants I : i\) G [0, 7 t/2 ], II : -0 e [7 t/2 ,7 t], III : ip G [7T,37 t/2 ] and IV : i/> G [37t/2,27t] 
then 


and moreover 


^ = n- 1 
= iz 1 


f dU 2 2 r ^ , _. 

mi-ii+m+iv a W 71 1 + r 

f dH J: 2 2 r „ .. 

/ dD^ = ^rT^ Ccos(0+5) ’ 

'-i+ii+m-iv a r 771 + r 


(35a) 


v* = n 1 


Qj, = 1Z 


f dn 2 2 r 

/_ I _„ +ra+lv #^ = ^rv^ Csm(0 - i) 

'—1+11+111—IV a V 7T i + r 


(35b) 


Note that all four symmetries have an upper bound 2/n. For known C, one may extract r 
and (j) ± 5 from these two sets of asymmetries. 


C. Statistics and Sensitivity Enhancements 

Before proceeding to numerical computation of C, let us pause to consider the statistical 
confidence in the extraction of r and (j). We focus on their extraction from the asymmetries 
(35). These asymmetries are expectation values of a random variable defined to take the 
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values ±1 on the quadrants I, II, III and IV, as specified in eqs. (35). For a sample of A 1 
events, the corresponding error ax = \/(l — X 2 )/N ~ 1 /\/A, for A" = and 

fl^,. The statistical conhdence at which one rejects the SM values Asm - thus measuring 
NP effects - is then characterized by the chi-square statistic (A — Asm) 2 /Ay- 

As shown in Ref. [13] and below, special ‘sensitivity parameter’ kinematic cuts may en¬ 
hance C on the resulting remaining phase space. The construction of these cuts is motivated 
by the observation, from eqs. (23), that C is enhanced on those areas of phase space in 
which the BH interference term, ~ a v 2 a 2 i , is comparable or larger than terms in the to¬ 
tal BH square amplitude |ck| 2 = J2\ rs l a rs| J - For example, one may define the sensitivity 
parameter 

$ — l Q h2 C %l/l Q hll 

~ 2(1 — cos[|</> + — (f )-1 — 7r]) 

Requiring S > 1 produces an event-level kinematic cut that typically leads to Of 1) en¬ 
hancements of C on the remaining phase space, as will be verified below. Note that the 
(1 — cos[|0 + — (f)- 1 — 7r]) factor in eq. (36) implies that S ~ 0 for coplanar events, and mo- 
roever that a lower bound on S typically favors events with higher acoplanarity. One may 
also consider other sensitivity parameters, such as 

T = 2|a 12 a 2 il/M > (37) 

which is normalized such that T E [0,1]. 

Let us define Co (Ao) to be the relative oscillation amplitude (number of events) in the 
absence of S or T cuts, and write 


7+7- 


(7+ + 7-) 2 


7+0+7- 

1(7^2 72^2), 


(l + 7^)(l + 7 OP). (36) 


C c = C[S >S C ,T>T C ], N c = N[S >S C ,T> %]. 


(38) 


Compared to the S c = % — 0 case, the application of sensitivity parameter cuts scales the 
NP statistical conhdence by the factor 


£ = 


A 


— Asm \' 

Ay J 


S>S C ,T > 71 : 


<JX 


A-Asm 


Ac 

A 0 


(39) 


for A = Q,p, and fThat is, the enhancement C c /Cq achieved by the sensitivity 
parameter cuts competes with the corresponding increase in the statistical error, A 0 /A c , 
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since necessarily N c < N 0 . We shall see in the next section that there are choices of S c 
and % for which £ > 1. For the purpose of measuring NP effects, this is equivalent to an 
effective increase in the sample size N 0 i—> TiN 0 - an increase in the effective statistics. 


IV. SIMULATIONS 

A. Relative Oscillation Amplitude 

Extraction of the relative oscillation amplitude coefficient C is achieved numerically via 
Monte-Carlo (MC) generation of B —> (K* —y K + 7T~)( r y —>• e + e _ ) events according to 
the differential rate (31). We use the matrix element in eq. (A4), generated from linearly 
polarized photon amplitudes. Though it does not provide the same analytical insight as 
the matrix element generated from spinor-helicity methods (23), this matrix element is 
as numerically stable as the latter, and moreover, acts as a convenient cross-check of the 
analytic results in eqs. (23) and (34). 

For the numerical results shown in this paper, we use a private MC code written in 
C/Python. For simplicity, we apply the narrow K* limit (7), which fixes the K* to be 
on-shell. We assume the nominal Belle II parameters [16] (see also App. B) in order to 
determine the B rapidity distribution /#(?/), which is peaked at /Ty ~ 0.29. More details of 
the operation of this MC generator are included in Appendix C. We have further checked 
the numerical results with a second private MC, written in C++/Java, that makes use of 
the matrix element (23). In both codes, we discard the overall normalization of the matrix 
element - e.g. A(r) in eqs. (23) - since we are concerned only with the relative oscillation 
amplitude, ~ C. 

In order to account for limited angular resolution, we include hereafter cuts on the lepton 
polar angles, 9±, and opening angle, 9u, defined in Fig. 1. We will consider a uniform polar 
cut 

Qa,± > @c j (40) 

for various values of 9 C . In particular, we consider two benchmark cases 9 C = 10 -6 
and 5 x 10~ 4 . The former captures almost all conversion leptons in the B —> ( K* —>■ 
K + 7r _ )(7 e + e _ ) process for semirelativistic B' s, while the latter might be plausi¬ 

bly achievable in the near- to mid-term future. To extract C, we fit eq. (34) to the 
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FIG. 2. Left: The fit value for C with the Tier error band as a function of the polar angle cuts 
6u,± > 0 C (see Fig. 1). The peak value of C approximately coincides with the peak of the 9± 
marginal distribution (see the left panel of Fig. 6). Right: Normalized differential distribution 
dIZ/dip for four different (r, </> + 5) couplets and 9 C = 10 -6 . Also shown are theory predictions 
(gray) for the input values of (r, cj> + 6) and the extracted value C[9 c = 10~ 6 ] in eq. (41). 

dlZ/dip histograms for various choices of r and (j) + 5, including the couplets {r, </> + <5} = 
{(0.1,0), (0.2, 7 t/4), (0.5, 7t/3), (1.0, 7t/ 2)}. In Fig. 2 we show the extracted C as a function 
of the 9 C cuts. The maximal relative oscillation amplitude one can expect is of 0( 20%), and 
the benchmark extracted C values are 

C[6 c = 10 -6 ] = 0.173 ± 0.001, and C[9 c = 5 x 10“ 4 ] = 0.150 ± 0.003 , (41) 

where the errors are purely statistical in origin. We also show in Fig. 2 typical dlZ/dip 
histograms for various choices of r and (fr+5. The expected shifted and amplitude-modulated 
cosine can be clearly seen. Applying the extracted value for C, these histograms are in 
excellent agreement with the theory predictions (34). 


B. Statistics Enhancements 

Incorporating the S and T kinematic cuts, we show in Fig. 3 the absolute enhancement 
of C as a function of the net cut efficiency, e, for the two benchmark polar angle cuts. The 
net cut efficiency is defined hereafter to be the fraction of events kept after application of 
both kinematic and polar cuts. The pure S cuts - that is, C(S > S c , T > 0) - provide the 
larger enhancement at high cut efficiencies. At lower efficiencies the pure T cut provides the 
larger enhancement. The C and e dependence on S c may be read off from Fig. 4. 
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FIG. 3. Upper panels: The coefficient C as a function of the cut efficiency for the S (blue) and 
T (gold) kinematic cuts, with polar cuts 0 C = 1CT 6 (left) and 6 C = 5 x 10 -4 (right). The colored 
regions depict the ±2cr statistical error bands. The equivalent effective statistics curve £ = 1, i.e. 
C = C 0 /^NjN~ 0 , is also shown (gray). Lower panels: Statistics enhancement £ as a function of 
the cut efficiency. The maxima correspond to the optimum cuts S > S° pt and T > 7^° pt (colored 
dots in all panels). 

Comparing to the equivalent effective statistics curve E = 1 - i.e. C = C 0 /\/N/N 0 - 
we see in Fig. 3 that the C dependence initially rises faster than Cq/^/N/Nq. This means 
that for low values of S c and %, the statistics on the S and T cut phase space is enhanced 
compared to the full data set. The explicit E dependence is also shown in Fig. 3. In 
particular, one sees that for 9 C = 5 x ICC 4 , the statistical enhancement E is optimized at 
5° pt = 1.1 ( T c opt = 0.62) corresponding to E = 3.6 (E = 2.5). 

To demonstrate the efficacy of these statistics enhancements, in Table I we compute 
the extracted values for r and <p + S, along with their statistical errors, with and without 
the optimal 5° pt and l~° pt cuts. We use an MC sample generated from an input couplet 
(r, 0 + 5) = (0.2, 7 t/ 4), containing a total of iV = 10 4 events. This roughly corresponds 
to 50 ab -1 of data - a benchmark luminosity at Belle II [16] - and a percent level photon 
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FIG. 4. The enhancement in C as a function of S. The secondary y-axis shows the corresponding 
cut efficiency. The colored regions depict the ±2 a statistical error bands. 

conversion rate with ideal acceptance efficiency. The polar cut is 9 C = 5 x ICC 4 . The 
extracted values are obtained by two different methods, first from a fit to the differential 
rate eq. (34) and second from the quadrant asymmetries (35). 


Method 

S,T> 0 

5 > 5 c opt 

T > T? pt 


r 

(0 ± (5)/7T 

r 

(0±<5)/7T 

r 

(0±<5)/7T 

dlZ/dilj 

0.194 ±0.017 

0.255 ±0.026 

0.194 ±0.010 

0.263 ±0.015 

0.190 ±0.013 

0.247 ±0.020 


0.217 ±0.005 

0.216 ±0.157 

0.203 ±0.003 

0.248 ±0.088 

0.192 ±0.003 

0.232 ±0.112 


TABLE I. Extracted values of r and (0 + 5) from an MC sample with input values (0.2, 7t/4). 

The statistical errors for the optimized kinematic S and T cuts are a 0(1) factor smaller 
than for the full data set, as expected from the above numerical analysis: Application of 
these sensitivity parameter cuts improves the statistical power of the r and 0 extraction 
from B —> (K* —>■ K + 'ir~)('f —>■ e + e _ ) events. Moreover, even without these enhancements, 
both r and 0 are extracted with sufficient numerical precision to probe NP effects at the 
r ~ 1 or 0 ~ 1 level. 

V. CONCLUSIONS 

In this paper we have presented the hclicity amplitudes and differential rate for the 
B —>■ K *7 —> K + 7 T “7 process, in which the photon undergoes subsequent nuclear conversion 
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to a lepton pair. Interference between the intermediate, on-shell photon polarizations in 
the coherent B — y (K* — > K + n~)(^ — > e + e _ ) process produces oscillations in the angular 
kinematic observable ip. Measuring the amplitude and phase of these oscillations - or equiv¬ 
alently two quadrant-type asymmetries and - permits extraction in principle of the 
polarization ratio, r, and the relative weak phase, 0, of the right- and left-handed B —y K *7 
amplitudes. In this manner, SM expectations for both r and 0 may be tested. 

We have employed private Monte Carlo simulations to compute the ip distribution and 
asymmetries as a function of r, <p and kinematic cuts. In particular, kinematic cuts with 
respect to the sensitivity parameters S and T may sufficiently amplify these oscillations, 
such that the overall statistical power of the r and <p> measurement is increased by an 0(1) 
factor. 

Implementing this approach using converted photons will be experimentally challenging, 
not least because of the high angular resolution required to reconstruct the conversion lepton 
kinematics. Moreover, a detector whose thickness is on the order of one radiation length 
or less is required to avoid multiple leptonic rescatterings, that otherwise smear the lepton 
kinematics. Nonetheless, the theoretically clean nature of the r- and 0-sensitive observables 
presented in this work may perhaps encourage the use of this technique in a future dedicated 
detector element. 
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Appendix A: The B —> Kire + e squared matrix element by polarization decomposi¬ 
tion 

The effective Lagrangian in eq. (2) gives the following Feynman rule for the BK *7 vertex: 

< ■7 (*0\m 

= i [g\\ {Am 2 BK *g^ - 2 k u p p ) + 2 g±e flUf)C7 k p p rT ] (Al) 

K*(p) K ,v 

where A m 2 BK * = m 2 B — m 2 K *. The amplitude for B — > 'y x Kn with A being the photon hclicity 
is then given by 


■A/a = 9k* { 9\\ [A m\ K *e\ ■ (p K - p n ) - 2 e* x ■ p [k ■ (p K - p n )}] 


- 2g±e tiUpcr ( p K - p 7r )^ £ u *k p p a 


(A2) 


The BH squared amplitude in the nuclear rest frame for a linearly polarized photon in the 
+z direction with polarizations A = { 1 , 2 } is 


bh av ~ ^ / gXW B q2 + ' p ~ + k • P+ ^ —if Ep+P ~ + Ep - pX A 

q A I (k ■ p~)(k ■ p + ) \ k ■ p_ k ■ p + J 


Ep+P- E p _p\ 
k ■ p_ k ■ p + 

(A3) 

where terms of (D(q 2 /E±) were dropped (see Appendix C in [13] for details). The squared 
amplitude is then given by 


\M\ 2 = Y AaAa ! /BH aa \ 

A,A'e{l,2} 


(A4) 


A numerical comparison between the above expression and eq. (23) shows excellent agree¬ 
ment over the entire phase space (sampled uniformly). 

Appendix B: B rapidity distribution 

Consider an e + e _ — y T — y BB factory, and let 0 denote the polar angle of the B 1 s with 
respect to the electron beamline in the center of mass frame. The amplitude for production 
•ATprod ~ sin©, and so the probability distribution 


Pn {cos 0) = - [l - cos 2 0] . 


(Bl) 
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Here and in the following we neglect effects of lab frame angular acceptance cuts, which may 
non-trivially restrict the domain of both rj and 0. 

In the center of mass frame - the rest frame of the T - each B has energy E* = mx/2: 
hereafter the * superscript denotes center-of-mass frame quantities. The corresponding ra¬ 
pidity, which we choose to be positive by convention on the branch 0 e [0, n\, 

rf = cosh _1 (mx/2m B ), (B2) 

and the B speed in this frame f3* = tanh(?/*). 

At H-factories the lab frame electron and position beam energies, E±, are asymmetric, 
but are chosen such that the T is on shell, i.e. 4£ + £_ = m\. For example, at Belle II the 
beams are planned to be £ + = 7 GeV and £_ = 4 GeV. The boost rapidity of the center of 
mass frame with respect to the lab frame is correspondingly 

r/ T = cosh -1 [(£ + + £_)/m T ] • (B3) 


The rapidity of the B in the lab frame may now be written as a function of cos 0, viz. 


r/(cos 0) = cosh 1 | cosh rf [ cosh r/x + ft* cos 0 sinh r/x] j , (B4) 


and its pdf, by definition 


/b(»?) = / dcos 0 po(cos 0)(i[r/ — r/(cos 0)]. 


(B5) 




Under a change of variables ( = rj( cos@), one finds 


Ib{v) = 


cosh rj* J\ vr - v *\ 


VT+V * d( ^ ~ ^ sinll C 


sinh rjr 


1 - 


/ cosh ( — cosh rf cosh r]r \ ' 
\ sinh if sinh r/x / 


3 sinh r/ 


4 sinh rf sinh r/x 


1 - 


cosh Tj — cosh rf* cosh r/x 
sinh rj* sinh r/x 


2n 


\vr —v‘\<V< r l"+Vr, 

(B6) 


and zero otherwise. Note that jf itself has zeroes at each end of its non-trivial domain, i.e. 
at r/ = rf + r/x and |r/x — rf\. The boost at the pdf peak is 7 pea k = (£+ + £-)/(2nr). E.g. 
for the Belle II parameters, the peak /37 pea k = sinh cosh _1 (7 pea k) = 0.29. This is the boost 
of B’s emitted at 0 = 7 t/ 2, and matches the quoted B design boost at Belle II [16]. 
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FIG. 5. Left: The cumulative distribution function (CDF) of the opening angle between the 
leptons. Right: The normalized distribution of the photon energy in units of the B mass for two 
different 9^ cuts. 

Appendix C: The Monte Carlo event generator 

This appendix describes in more detail the MC event generator written in C and Python. 
The B —y (K* —y Kn)^/ phase space is generated as follows. The B rapidity is sampled from 
the PDF given in eq. (B6) while the photon polar angle # 7 and the K polar and azimuthal 
angles are generated uniformly in the appropriate frame. On the other hand, since in BH 
photon conversion the leptons are produced with preferentially small angles with respect to 
the photon direction, the lepton polar angles are generated uniformly on a log scale. This is 
implemented via the transformed variables t± = log 10 9± where t± are uniformly distributed 
and with t± G [—5, —1]. Moreover, the azimuthal angle separation between the leptons (5<fi) 
is peaked around n and so, to improve the efficiency of the generator, 8<p is sampled from a 
Cauchy distribution. All other Bff variables are generated uniformly. 

The weight associated with each event is proportional to the matrix element (A4). The 
events are unweighted using the standard procedure. That is, the weights are normalized 
to the largest weight and the event is kept if its normalized weight is larger than a random 
number on [0,1]. Of course, this procedure assumes that the phase space was sufficiently 
sampled such that the largest weight found is close to the global maximum. 

Using this procedure, we generate MC samples for many choices of (r, </> + 5) couplets 
with 500k events per sample. Some representative distributions from the (0.1,0) sample 
are shown in Figs 5 and 6. In particular, the left panel in Fig. 5 shows the cumulative 
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FIG. 6. Left: The normalized polar angle distribution of the positron for two different values of 
the opening angle cut 9^. Right: the positron energy as a fraction of the photon energy for two 
values of 6y>- The distribution exhibits the expected behavior for BH conversion. It is symmetric 
about 1 /2 and prefers that one lepton carry a larger fraction of the photon energy. 


distribution function for the opening angle between the leptons 9u while the right panel 
shows the distribution of photons energies. Figure 6 shows the polar angle and fractional 
energy distribution of the positron. 
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